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TECHNICAL NOTE 3956 


LEFT AND MOMENT RESPONSES TO PENETRATION OP SHARP-EDGED 
TRAVELING GUSTS, WITH APPLICATION TO 
PENETRATION OP WEAK BLAST WAVES 
By Joseph A. Drischler and Franklin W. Diederich 

SUMMARY 


The lift and moment responses to penetration of sharp-edged 
traveling gusts are calculated for wings in incompressible and supersonic 
two-dimensional flow, for wide delta and rectangular wings in supersonic 
flow, and for very narrow delta wings. By using the two-dimensional 
Indicia!-- ] 1 f t functions, some calculations of normal-acceleration response 
are made for two mass ratios. 

The results of these calculations indicate that the forward Bpeed 
of the gusts has a large effect on the lift- and moment-response func- 
tions. For incompressible flow, peaks exist in the early portion of the 
lift response, which may be much larger than the steady-state value. 

Some peaks also occur in the lift-response functions for supersonic 
speeds but are much less pronounced and exceed the steady-state value in 
only a few Instances. 

Calculations have also been made of the normal-acceleration response 
to sharp-edged traveling gusts and indicate that this response tends 
to follow the lift response very closely in the first few instants of 
penetration; thuB, the large peaks which exist in the lift response at 
subsonic speeds are duplicated in the acceleration response. 

The relation between gusts traveling at supersonic speeds and blast 
waves is indicated, and the manner in which the calculated lift and 
moment responses can be used in a linearized approach to the blast-load 
problem is outlined. 


INTRODUCTION 


The growth of the lift and moment on a wing entering a stationary 
sharp-edged gust has been the subject of numerous investigations since 
it was first calculated for incompressible two-dimensional flow in 
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references 1 and 2 . However, very little work appears to have been done 
on the subject of lift and moment response to traveling sharp-edged gusts, 
the only published results being those presented for incompressible two- 
dimensional flow in reference 3. 

The reason for this shortage of information 1 b probably due to the fact 
that gusts (in the literal sense) traveling at high speeds are not likely 
to be encountered in practice. However, the lifts and moments due to 
penetration of traveling gusts can be used to calculate the loads »nr) 
motions of an airplane or missile dropped from another airplane, of an 
airplane crossing the wake of or flying past another airplane, of a 
helicopter blade traversing the wake produced by itself end the other 
blades of the rotor (a problem which furnished the motivation for ref. 3), 
of an airplane flying through a sonic "boom," and of an airplane encoun- 
tering or being overtaken by a blast wave. 

In the first section of this paper the unsteady lift and moment 
pursuant to penetration of a traveling sharp-edged gust are calcu- 
lated by linear theory for incompressible two-dimensional flow (by a 
method different from the one employed in ref. 3), for two-dimensional 
supersonic flow, for delta wings with supersonic leading edgeB, for 
recta n gu l ar wings in supersonic flow, and for very narrow delta wings 
in incompressible and compressible flow. The unsteady- lift functions 
presented here may be considered to be generalized unsteady- lift functions 
in the sense that they Include as special cases the two previously calcu- 
lated unsteady-lift functions: namely, the gust-penetration (Kiissner) 

function, which corresponds to zero gust speed, and the unsteady- lift 
function due to airplane motion (the Wagner ^function, first calculated 
in ref. k), which corresponds to infinite gust speed. 

By using the -unsteady-lift functions calculated for two-dimensional 
incompressible and supersonic conpressible flow, the normal-acceleration 
response is calculated in the second section of this paper for airp lane s 
with two different mass ratios and for various values of the speed of 
gust propagation in order to indicate the effect of this speed on the 
accelerations . 

The relation between traveling gusts and blast waves is indicated 
in the third section of this paper. The results presented herein for 
supersonic speeds are thus pertinent to the penetration of blast waves, 
although the results calculated by linear' theory are usable only if the 
blast waves are weak, and although for blasts the unsteady-lift functions 
calculated in this manner represent only the effects of the change in flow 
direction' induced "at the wing. The effects of the finite pressure Jump, 
which are of first order, are not taken into account entirely w- nfl may 
require separate treatment; and neither are the effects of the discon- 
tinuities in density and temperature taken into account. However, these 
effects are of higher order and, hence, negligible for a weak blast wave. 
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General relations "between the lift or moment responses to sharp- 
edged gusts (traveling or stationary) and, the responses to indicia! flap 
deflections are given in the appendix. These relations have "been used 
to calculate the unsteady- lift and unsteady-moment functions given herein 
for incompressible two-dimensional flow and to check the results obtained 
for some other conditions. 


SYMBOIS 


A 

a 

b(x* ) 
C(k) 

C m 

CWCa) 

°P 

c 

C Z 

c, 

4 a 

c m 

d 

k 

K 

k 


aspect ratio 

speed of sound of undisturbed stream, 
local span 
Theodorsen function 
lift coefficient 

slope of lift curve 

pitching -moment coefficient (steady-state value) 
instantaneous value of pitching -moment coefficient 
pressure coefficient 

chord, root chord in case of delta wing 
section lift coefficient 

section slope of lift curve 

section pitching-moment coefficient 

distance behind wave front 

local effective sinking speed of wing 

alleviation factor for sharp-edged gust 

reduced frequency, — 
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k(s) 

unsteady-lift function for 
edged gusts 

k 8 (s,|) 

unsteady- lift function for 

k^s) 

unsteady-lift function for 

^(s) 

unsteady-lift function for 

Z(x') 

lift per unit length along 

M 

Mach number, v/a 

«e 

equivalent Mach number, M 


penetration of traveling sharp- 

indicial flap deflection 
indicial change of angle of attack 
penetration of stationary gust 
chord _ _ 

+ M g 


“g 

m 

“a 

Pb 


Mach number corresponding to gust propagation speed 

slope of leadin g edge, b(x')/2x' 
mass of airplane 

differences in pressures on upper and lower wing surface at 
a given point 

overpressure in blast 


<1 


dynamic pressure, 


Ipv 2 


S wing area 

s distance traveled, semi chords or root semichords (unless 

specified otherwise) 

t transformed time (distance traveled by sound wave), at' 

t’ physical time 


V forward velocity of wing 


V g propagation velocity of gust, positive when moving in 

s direction opposite to that of wing 


W 


vertical velocity of gust 
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. 

1 


x 

X* 

Z 

2b 

“o 


transformed ordinate corresponding to x 1 , ^x' - Vt*) 

ordinate along flight path, measured rearward from leading 
edge or apex of wing 


acceleration of airfoil due to traveling sharp-edged gust 

qsc, 

■Uv w 

static acceleration, — 

% V 


angle of attach of wing prior to gust entry, radians 


3 -v 1 1 - m 2 ! 

8(b) 

A 


unit Impulse function or Dirac delta function 
V 


speed ratio. 


mass ratio. 


V + V, 


e 


H. 

pScC L a 


p air density 

0 perturbation velocity potential 

as circular frequency 

_7(s) unit jump function 

Subscripts : 


b blast 

t time 


Dots indicate differentiation with respect to t ime . 
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CALCULATION OF UNSTEADY- LIFT AND UNSTEADY-MOMENT FUNCTIONS 
FOR TRAVELING SHARP -EDGED GUSTS 


The unsteady-lift and unsteady-moment functions for indicial angle- 
of -attack chang e and for penetration of a sharp-edged stationary gust 
have been calculated for a large variety of conditions. (See refs. 1 
to 11, for instance . ) This information is not sufficient in itself to 
calculate the unsteady-lift B.nd unsteady-moment functions for traveling 
sharp-edged gusts.. Therefore, in this section most of the functions are 
calculated directly by solving the given boundary- value problem by use 
of the te chni ques employed previously to solve the other problems. 

The desired functions can be calculated if additional information 
is available for the indicial-response functions, such as the Indicial 
responses to flap deflections. This approach is outlined in the appendix 
and will be used in this section to calculate the deBired responses for 
incompressible two-dimensional flow. 


The Boundary-Value Problem 

The linearized partial differential equation for the velocity 
potential of a wing traveling at a constant speed V in a compressible 
inviscid medium is the wave equation and can be written in coordinates 
fixed in the wing surface as 

(l - m2 )0x'x' + 0y'y' + 0z’z’ = ^(^x't’ + 0t’t') C 1 ) 

BT 

and for the lifting case the pressure difference £p between the upper 
lower surface is given in terms of the potential on the surface by 

f - * 0t') < 2 > 

For the sharp-edged gust-penetration problemjthe primary boundary con- 
dition is that the downwash 0 Z , be zero on the part of the wing which 

has not yet been reached by the gust front, and be equal to W on that 
part of the wing which has been passed overly the gust front; that is, 
on the wing surface. 



MCA TS 3956 


7 


0 B * “0 (x> > (V + V g )t') 

0 Z> -W (x* < (V + Vg)t*) 

Other boundary conditions axe that the pressures on the plane of the 
wing he zero (except on the wing) and that the pressure difference at 
subsonic trailing edges he zero. 

For supersonic flow the following coordinate transformation often 
facilitates the solution of the problem: 



x - x' - Vt* 

y - y* 

z = z* 
t = at 1 

Under this transformation, equation (l) becomes the normalized wave 
equation in stationary coordinates 

0xx + $yy + fzz " 0tt (5) 

equation (2) becomes 

< 6 > 

ani equation ( 3 ) becomes, for the transformed wing surface. 



0z “ 0 ( x > V) 

0 Z - V (x < Mgt) 


(7a) 


for gusts traversing the wing from the leading edge to the trailing edge. 
Similarly, 
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0 Z =■ 0 (x < c + 

0 Z “ W (x > c + Mgt *) 

for gusts overtaking the wing. 


(7b) 


Definitions of Unsteady-Lift and Unsteady-Moment Functions 

The section lift and moment coefficients can he obtained by inte- 
gration of the pressure distribution along a chord, namely. 





On a wing of finite section the loading on a spanwiBe strip can be 
obtained by integrating over y, namely. 


q. 


I. 


b(x ' ,/2 ^a y . 

-b(x')/2 % 


(9) 


and, hence, the lift and moment coefficients 


can be obtained from 



( 10 ) 
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The unsteady- lift functions for penetration of traveling sharp- 
edged. gusts will he nor mal ized, to their steady-state values and, hence, 
he defined for two-dimensional flow as 


k(s) = 


c 7 (s) 




and for three-dimensional flow as 


k(s) s 


Cl/ 8 ) 

C H 
%.V 


where both c^Cs) and Cj^s) pertain to vertical gusts of intensity W, 

and s is the distance traveled in semichords (for two-dimensional and 
rectangular wings) or root semi chords (for delta wingB). 

For the unsteady-moment functions no common designation exists. 
Therefore, unsteady-moment coefficients will he presented directly for 
unit W/V and will not he normalized with respect to the steady-state 
values. 

When the gust propagation speed is zero, the function k(s ) reduces 
to the commonly used gust -penetration function, which is usually desig- 
nated hy Is^s) and referred to as the KluBsner function. Thus, 


¥b)| v ^ 

8 

Similarly, when the gust propagation speed is infinite, the angle of 
attack over the entire wing is changed instantaneously; therefore, the 
unsteady- lift function becomes that for indicial change of angle of 
attack, which is usually designated hy k]_(s) and is referred to as the 
Wagner function. Thus , 


h(s)| y =oo = k-^s) 
6 


Between these extremes, when the gust propagation speed is positive, the 
gust approaches the. wing; when it is negative hut less in magnitude than 
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the speed of the wing, the wing overtakes the gust. When the gust speed 
is negative hut greater in magnitude than that of the wing, the gust over- 
takes the wing; that Is, it approaches the wing from the rear. All these 
possibilities will be considered in this paper. 


Incompressible Two-Dimensional Flow 

For incompres sible two-dimensional flow, equation (l) becomes the 
two-dimensional Laplace equation, and its direct solution for the uns teady 
boundary condition represented by equation ( 3 ) constitutes a dif ficult 
problem. This problem has been solved in reference 3 by an indir ect 
approach, which consists in using the known results for the lift response 
due to an infinite train of traveling sinusoidal gusts and in obtaining 
the response for sharp-edged gusts by means' of the well-known super- 
position integral given in reference 5 * However, in order to effect 
this transformation, the results for sinusoidal gusts have to be expanded 
in a series; thus, the results for the sharp-edged traveling gust con-, 
tain a certain degree of approx imat ion. 

The same function had been obtained by an altogether different 
method in connection with the present paper before reference 3 became 
available. This method is based on the approach outlined in the appendix 
and consists in relating the lift or moment response to penetration of a 
stationary or traveling sharp-edged gust to'indicial lifts and moments 
due to flap deflection by means of superposition integrals. Its appli- 
cation to the case of incompressible two -dimensional flow is outlined 
in the following paragraphs. Although this, approach. also contains an 
approximation, a comparison of its results with those of reference 3 
should furnish an indication of the validity of the two approximations. 

Equations ( Al) and ( A2) of the appendix serve to express the desired 
unsteady-lift function k(s) for traveling gusts in terms of the unsteady- 
lift function for indiclal flap deflection kg(s,|) for the wing with 

the given plan form and Mach number; for incompressible two-dimensional 
flow this function kg 1 b given by equation (aB). Thus, upon substi- 
tuting this function into equation (Al), the following expression is 
obtained for traveling gusts which pass over the wing from the leading 
edge toward the trailing edge: 






2* 


8 (s-Aff) 




dcr 


(ID 
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where k-^ c (s ) is the continuous part of the 1 ndl cial-llft function 
k]_( s ) , that is, the part of kj_(s) excluding the delta function for 

this case. (The notation in the upper limit of the integral specifies 
that the upper limit is the smaller of the two quantities 2 and s/X. ) 
The functions T^ and T^q are (see eq.. (A6)) 

Tj,. = (| - 1)\(2| - | 2 - cos'Hi - 

T 10 = \] 2* - % 2 + cos-^g - 1) 

Similarly, for gusts overtaking the wing (passing over the wing 
from the trailing edge toward the leading edge), equations (A2) and (a 8) 
yield 



k(s) 


= /o 






iki >c (s-|x|a) 


-dT. 


10 


|=2-ct 


1_ 
2 it 


8(s-|A|a) 


r-cffi], 

r*r 


g=»2 -a 


da 


(13) 


Hence, for gusts approaching the wing or receding from it at a 
speed "below that of the wing, equations (ll) and (12) yield 


k ( 8 ) " i J q ^ k l( B -^ B l)\fl^ ^1 + -^2 \I B ^ - ®j 

k ( B) =if 0 


(a ^ 2X)] 


>( 1 ^) 


(s £ 2X) 


J 
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and, similarly, for gusts overtaking the wing, equations (12) and (13J 
yield 


■ s/I * 


p a/ | a | 2-8, 

k(s)=£ ^(s-IXlsJJ—- ds 1 + 


jrA 


~ \js(2|?v| - s) 


(b < 2| >T|) 


(15) 


2 

k(s) = J k l(s~l ^l B i)v 2 — (s>2|x|) 


In the same manner, expressions for the unsteady-moment function 
Cj^s) for traveling gusts can he obtained from equations (Al) and (A2) 

by using the moment function for indicial flap deflection given in 
equation (All) and the function Tg defined in equation ( A9 ) » 


The expressions for the lift due to penetration of traveling gusts 
cannot be evaluated exactly because no explicit expression for k-^ c (s) 

is known. However, the following simple function was shown in reference 5 
to be a good approximation to the Wagner function: 



<■) 


** 1 - 


2 

4 + 8 


(16) 


With this expression the unsteady-lift functions k(s) can be evaluated 
readily. The results are given in table 1 and are plotted in figure 1 
for several values of A. (in fig. 1 the abscissa is interrupted at 
s *= 5 and different scales are used in order to exhibit the behavior of 
the response functions at both small and large values of b.J For A = 0 
(Vg = ± »),, the functions k(s) given by equations (l4) and (15), respec- 
tively, reduce to_the approximate Wagner function given in equation ( l6)j 
and for A = 1 (Vg =* 0), the expression for gusts approaching the wing 

reduces to an equivalent approximation to the Kussner function. 

The analytic expressions given here for k(s) can be compared 
directly with those given in reference 3 only for small and very large 
values of s, in which range the agreement J.s very good. In the inter- 
mediate range of s the numerical results presented in figure 1 agree 
with those estimated in reference 3 to within the width of the lines. 
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Therefore, In view of the entirely different nature of the approxi- 
mations made in the two methods, the excellent agreement of the results 
is gratifying. 


An examination of figure l(a) indicates that, as the propagation 
speed of the gust increases from 0 to 00 (x decreases from 1 to 0), 
the continuous part of the lift-response function k(s), which represents 
the circulatory contribution, increases gradually for an finite val ue s 
of s from the Kussner function to the Wagner function. At the same 
time the noncirculatory contribution increases also and more rap idly , 
giving rise to a hump in the response function with a maxi mum at approxi- 
mately the value of s which corresponds to the t ime it takes the gust 
to traverse one -half of the wing chord. This hump, which exists neither 
in the Kussner nor the Wagner function (as usually defined), may be very 
large, much larger than the steady- state value, in fact. Ab the forward 
speed of the gust approaches infinity, this hump approaches the function 


^ 8(s), which represents the noncirculatory contribution associated with 


the singular part of the Wagner function. 


When the gust recedes from the wing at speeds less than that of the 
wing (A > l), both the circulatory and the noncirculatory part of the 
response tend to decrease until, when the gust travels almost as fast 
as the wing, the response increases very slowly und thus takes a very 
long time to attain its steady-state value of unity. 

The lift-response functions for gusts overtaking the wing, that 
is, for Vg < -V and, hence, for negative values of A, are shown in 
figure 1(b). For very large gust speeds the response function again 
approaches the Wagner function. For smaller speeds both the circulatory 
and the noncirculatory part of the lift gradually decrease, until, when 
the gust overtakes the wing with a speed barely higher than that of the 
wing, the lift increases very slowly toward its steady-state value. As 
in .the case of gusts approaching the wing, the noncirculatory contri- 
bution again leads to a hump which may be very large. 

The circulatory part of the lift tends to be greater for gusts over- 
taking the wing t han for gusts approaching the wing or receding from it, 
for the same reason that tral ling-edge flaps are more effective in pro- 
ducing lift than leading-edge flaps, although the noncirculatory part 
is the same (for a given value of | A|). Thus, for a gust overtaking a 
wing the lift response is larger at all finite values of s than for a 
gust being overtaken by the wing at the same relative speed, namely, 

l v + v s!- 

The unsteady-moment functions c^s) which were calculated in a 
similar manner are given in table 2 and are shown in figure 2. As may be 
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seen from equations (A10) and (All) by the absence of terms proportional, 
to C(k) and k^^ 0 (s), the moment about the quarter- chord points con- 
tains no contribution due to circulatory effects. Hence, no approximation 
of the type indicated in equation ( 16 ) needs to be made for the moment 
function. The behavior of the functions c^s) shown in figure 2 is 
thus similar to the variation of the humps in k(s) discussed previously. 

As the forward propagation speed of the gust approaches infinity, 
the moment function approaches the function -jt5(s), which indicates 

that the lift represented by the contribution Jjj B(s) to k(s) for that 

case acts at the midchord. As the forward speed of the gust decreases, 
the initial center of pressure moves forward until it reaches the quarter- 
chord point of the wing when the forward speed of the gust is zero. As 
the gust recedes from the wing at velocities 0 > V g > -V, the initial 
center-of -pressure location is forward of the quarter chord, thus 
resulting in positive initial pitching moments. 


Supersonic Two-Dimensional Flow 

For supersonic two-dimensional flow the functions k-^( s ) and ^(s) 

have been calculated in references 6 and 7* The function k(s) and 
Cjj^s) can be obtained in the same manner. For this purpose equations ( 5 ), 
( 6 ), and ( 7 ) are convenient, because they serve to identify the actual 
boundary-value problem ( in the transformed coordinates ) with the boundary- 
value problem of a three-dimensional wing with supersonic edges in steady 
flow parallel to the t-axis with a Mach number of \j~2. The potential 
and pressure for this wing are then also the ones for the problem of 
interest here. 


The fictitious three-dimensional wings are shown in figure 3* The 
leading and trailing edges are indicated, and the leading edge also 
includes the part of the t-axis between 0 ..and c. The boundary con- 


ditions are that in the shaded region 




-H, 


elsewhere on the wing 


0 Z = 0, and ahead of the wing 0 = 0. For these boundary conditions the 
pressure coefficient can be obtained readily by a superposition of sources 
or by conical-flow methods, and these are given in table 3* The lift and 
moment coefficients as a function of time can be obtained from equa- 
tion ( 8 ), and the lift function can be reduced to k(s) by dividing by 

K k 

T ll^“ 


The results for k(s) and c^s) are given in tables 4 


and 5 and are shown for M = 2 in figures b and 5 for several values 
of Mg. For Mg = 0 end Mg = «> the function k(s) reduces to the 

functions k^(s) and given in reference as may be expected. 
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As nay be seen from figure 4(a), as the forward speed of the gusts 
increases, the rate of lift growth at small values of s increases and 
tends to approach infinity as the gust speed approaches infinity. For 
low gust speeds the increase is steady, but for high gust speeds 
^Mg ^ 1 for M = 2) the lift increases to the initial value of the k]_ 

function and dips down before it increases again toward the steady-state 
value of unity. The large peaks associated with the noncir culatory con- 
tribution in incompressible (and presumably subsonic) flow are not present 
for supersonic flow. 

For gusts receding from the wing the gust-response function increases 
steadily, and the rate of increase decreases with gust speed until, for 
gusts receding at almost the speed of the wing, the function increases 
very slowly. 

Similarly, as may be seen from figure 4(b), for gusts overtaking 
the wing, the lift response increases linearly (at a rate which increases 
with the speed of the gust) until it reaches a value intermediate between 
the initial value of the Wagner function and the steady-state value of 
unity, and then continues increasing at a much lower rate. 

Similar observations may be made for the moment functions shown in 
figure 5. 


Delta Wings With Supersonic Leading Edges 

For delta wings with supersonic leading ed ges, th at is, for delta 

wings with an aspect ratio A larger than k/W - 1 or a slope m 
of the lead ing edge (tangent of the apex half -angle) larger than 

1/yM 2 - 1, the functions k^( s ) and ^(s) have been given in refer- 
ences 8 and 9* Again, the functions k(s) and Gm( s ) can be obtained 
in a similar manner. 

The detailed solution of this three-dimensional-flow problem is much 
more difficult than that of the two-dimensional -flow problem discussed 
in the preceding section. However, if attention is confined to the chord- 
wise load-distribution function Z(x) and the total lift and pitching 
moment, the problem can be reduced to a two-dimensional one similar to 
the problem considered in the preceding section. This reduction can be 
effected (see ref. 8) by considering the integral of 0 over the local 
span, that is, 



<3y' 


<t> i= 
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for which. In the coordinates x,z,t, the partial differential, equation 
is equation (5)# with the y derivatives deleted. The loading func- 
tion Z(x) is then given by the equivalent of equation (6), n ame ly. 


itsl 

q. 



and the boundary condition is the equivalent of equation (7), namely, 
for the fictitious wing surface in transformed coordinates. 


$ 


z 


L 


m( x+Mt ) 
■m( x+Mt ) 


0z 


Thus, 


or 


<t z => 0 (x > Mgt) 

® z =» 3Wn(x + Mt) (x < Mgt) 


4 Z = 0 (x < c + Mgt) 

® z = 2Vfin(x + Mt) (x > c + Mgt) 


depending on whether Mg is greater than or_less than -M. 

The geometric characteristics of the fictitious wing are again 
those indicated in figure 3 (where c now refers to the root chord) and 
the bo undar y condition on f is again specified on the shaded part of _ 
the wing and is zero elsewhere; but it is not constant as it was before. 
The solution of the problem cannot, therefore, be obtained readily by 
conical-flow methods; but the source- superposition method can still be 
used conveniently, that is, 


<t(x,0,t) 


i II 

Mach 
fore cone 


*z (xi^lfol dt l 
]j(t - t ± ) 2 - (x - Xl ) 2 


The resulting expressions for the loading function are given in 
table 3, and the corresponding values of k(s) and C m (s) are given in 



L 


NACA TN 3956 


17 


tables 4 and 5 and are shown in figures 6 and 7 for M = 2 and several 
values of Mg. Again, the function k(s) reduces to the known func- 
tions k^s) and k2(s) (see ref. 9, for instance) for Mg = » and 
Mg = 0, respectively. The "behavior of these functions may he seen to 
he similar to that of the functions for supersonic two-dimensional flow. 
(See figs. 4 and 5*) 


Rectangular Wings in Supersonic Flow 

The pressure distribution in the region of a rectangular wing ahead 
of the Mach cones emanating from the wing tips is identical, to that of 
a two-dimensional wing. Although the pressure distribution within the 
Mach cones ca nn ot he calculated readily, the chordwise loading and the 
total lift and moment contributed by these regions can he determined 
readily by a technique similar to that used in the preceding section 
provided the Mach cone emanating from one wing tip does not cross the 

opposite side edge, that is, provided A is larger than iHFTL 
For such wings which will he referred to as wide rectangular wingB, the 
functions k^s) and kgCs) have been given in reference 10 and 
elsewhere . 

In reference 10 the unsteady- lift coefficient for a wide rectangular 
wing with a downwash distribution of the form 

0 z '(x , ,y , ,O,t«) = Vf(x' iZ(t*-t^) 

(where J_ represents the unit jump function and t^ may depend on x' ) 
is expressed as 


C L (t') = c L,o( t * ) " *Cl(*'> (17) 

where f s 1 ndlclal-lift coefficient for a two-dimensional, 

wing with the same downwash condition. 

The correction term which must be added to the two-dimensional value 
can be expressed as an integral of a chordwise -loading function 7 (x',t') 
and is, for a wing of unit chord, 


ac L ( t 1 ) 


k- 

a(m 2 - 1 ) 



(18a) 


An analogous correction for. the moment coefficient is 
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^M 2 


r f y(x' ,t' )x' dx' 

- l) J o 


(l8b) 


where 


2C ^ 

7 (x’,t') = f eu 
J 0 

— — 
x 1 - i 

f(6)dg ( 19 ) 

V(t’ - 

where. In turn, the indlcial-response 

function 6 l-j_( t) Is given hy 

a^T) =» 0 

(t < 0) 


- 1 

^0 = t = 


m l 

= T 1 

( m ! - T - ^ 

I 

= 0 

( T = “l) 

and 



_ _ M 

m l " “ 

- 1 
M 


M + 1 
n-, = ■ ■ ■ 

1 M 


Now, for a wide rectangular wing penetrating a traveling sharp- 
edged gust of intensity W and speed V g , 

f(S> - 

, a 

V 

(20) 

^ 's(r 

(•%>-“) 

(21a) 

and 



to ^ ) “ -r" £ 

a|Me| 

(Mg < -M) 

(21b) 



NACA 031 5956 


19 


where 


Mg = M + Mg 

Substitution of equations (20) and (21a) into equation (19) yields 



or, changing the variable of integration to 


r = 


= ^g. x' - 6 


M Mgt - (■ 


gives 


7(x, ’ t) = f I 0 


xl 

Mt 


*i(t) 




dr 




(22) 


and, similarly. 


r(x ’’ t) -f v( x ' 


MgX' 


- 1 


- M /o 


M(Mgt+l) & 1 (t) 


( x - *)’ 


dr (Mg < -M) (23) 


The geometric characteristics of this problem are indicated in 
figure 8. 'Phe integrations required in equations (22) and (23) take 
place along lines Mb = Constant from g = x’ (where t => 0) to {• = 

^wtoere t = ^ or t => 7 depending on whether Mg £ -MJ. 


0 
The 


resulting expressions for 7 are given in table 6. These expressions 
have been integrated to obtain the correction for the lift and moment 
coefficients. (See eqs. (l8).) The results are given in tables b and 5. 
The unsteady- lift function k(s), which is obtained by normalizing the 
corrected value of the two-dimensional unsteady-lift function (see 
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eq. (17) ), is shown in figures 9 to 12 for JJ = 2 and several values 
of Mg and A. Similarly, the unsteady-moment coefficients are given... 
in figures 13 to 1 6 for several values of Mg and A. Again, the 
unsteady- lift function k(s) reduces to the known functions k^Cs) and 
k 2 (s) for this case, which are given in. reference 10. 

The behavior of the lift-response and moment -response functions is 
similar to that for supersonic two-dimensional flow, except for the 
slight differences associated with the' f act !that for f ini te aspect ratios 
the initial, value of the k-j_ function can be greater than the steady- 
state value. 


Very Narrow Delta Wings 

For very narrow delta wings and bodies, derivatives with respect to 
x' or x tend to be small compared with the others in the equations of 
motion (eqs. (l) and (5))« For incompressible flow, equation (l) there- 
fore becomes the two-dimensional Laplace equation in the Y'Z' -plane anfl 
can be solved readily for both steady and unsteady boundary conditions. 

If only the chordwise-loading function and the total lift and moment 
are of interest, the solution of the boundary- value ' problem can be 
avoided, because the chordwise loading can be obtained directly from 
apparent-mass considerations . This approach^vas used in reference .11 
to calculate the functions k^( s ) and ^(s). 

For compressible flow, equation ( 5 ) becomes the two-dimensional- 
wave equation with boundary conditions in the case of interest here which 
can be construed as pertaining to a three-dimensional wing in steady 
supersonic flow, as was done for the functions k^( s ) and k 2 (s) in 
reference 12. 

Incompressible flow .- The expression for. the chordwise loading given 
in reference 11 is 


l(x«) 



where b(x)/2 is, for a delta wing, equal to mx', m be ing the tangent 
of the semiapex angle of the wing. For a wing flying at speed V and 
penetrating a traveling sharp-edged gust approaching from the leading 
edge with speed V , the effective sinking speed £ is given by 

o 



(2k) 
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where V e s V + V„ and J designates the unit jump function, as 
before. The total indiclal lift and moment about the two-thirds root- 
chord point of the wing are thus 


k(s) 


|L[l -2 (s-2M] + 1 (=-a) 


GmCs) 

C L (“) 


24x^ 


(3 - X) 


sf_ 

6x3 


1 -2 (s-2X) 


When X is zero or infinite, the following two limiting caseB are 
obtained: 


k(s) ix=i = T (s “ 2) 

= 1 (s>2) 


and 


k < s )U=o = f 6(b> + 1 (e> 

where 5(s) is the Dirac delta function. These results are identical 
with the functions k-^s) and kg(s) presented In reference 11 for 
this case. 

Similarly, for gusts overtaking the wing. 


hCxSt' ) = W 2 + 


1 . c ~ x ' 


Thus, the lift and moment functions become 


k(8) ‘ IxT l 1 ‘ Hxt) 2 & - 1 (s - 2lM) ] + 1 


Cm( s ) 

C L (*>) 


[2 - (1 + ^(3 -X) 

1 -2 (s— 2 I x!) 

L V 2 X/ J 

— J 
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(’(a) 

where — is again referred to the two-thirds root-chord point. For 

C L (“) 

the special case for which V g is infinite, this function k(s) reduces 

case 


to the function k-j_(s) given~in reference 11. When V g = -V (the 

of a gust traveling just behind the trailing edge), k(s) = 0; however, 
if = -11m (V Q + e) (the case of a gust traveling at the same velocity 


V fl = -11m 
" e-»0 


as the airfoil and with its leading edge just ahead of the trailing edge 
of the wing), then 


k(s) = _2( s) 

The functions k( s ) and C m ( s ) for gusts approaching the wing 
from either direction are shown in figures 17 and 18 for several values 
of 


The behavior of the response functions is similar to that of the ” 

functions for incompressible two-dimensional flow. Again, large peaks 1 

associated with the noncirculatory contribution to the lift exist at 
low values of s for gusts approaching the -wing or overtaking the wing; 
and, again, the value of the response function at all finite values of ... ,e 

s tends to be much larger for gusts overtaking the wing than for gustB 
approaching it or receding from it at the same relative speed. 

* 

Compressible flow .- For compressible flow the lift and moment 
responses can be obtained by the method used in reference 12. The 
resulting expressions are given in tables J and 8. They involve a func- 
tion f ( t] ) which is defined and tabulated in reference 8. Inasmuch as 
no simple explicit expressions for f( rj) are available, the integrations 
required to obtain the functions of Interest .here would have to be per- - - 
formed numerically. No such calculations have been made. However, the 
unsteady- lift functions given here can be shown to reduce to the func- ... IT 

tions k^(s) and ^(s) given in reference 12. ' 


NORMAL- ACCELERATION RESPONSE QF AN AIRPLANE 
TO TRAVELING GUSTS 


The purpose of this section is to indicate the effect that changes 
in the response function which result from changes in the gust speed 
may have on the acceleration response of an airplane to sharp-edged gusts. 
For this purpose the airplane will be assumed to be free to move in only 
one degree of freedom, vertical translation.' 7 '; Ignoring the pitching 
degree of freedom implies that not only pitching motions and their effect 
on the normal acceleration but also the unsteady moment which results 
from the fact: that the gust strikes the wing; and tail at different 


■ ftfl I 1 ■ . • : I lbs Jk I I 1 l.l Hli I li tl 
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instants of t ime are disregarded. Nonetheless , this assumption is made 
very often in gust-load studies, particularly in primarily qualitative 
studies, nnfl will he used here for that purpose. For airplanes with a 
Vi-i gh mass ratio, with a relatively high pitching moment of inertia, or 
with a well-damped short-period mode, or for literally sharp-edged gusts 
(such as those represented by a blast), this approach is even capable of 
yielding useful quantitative information. 

In terms of the unsteady- lift functions k(s) and k-j_(s) due to 

traveling-gust penetration arid plunging, respectively, the equation of 
motion for this case can he written in dimensionless form elb 

K(s) - k(s) ~ ^ f Q k l( s - s l)K(s 1 )ds 1 (25) 

where K(s) is an alleviation factor referred to the reference "static" 
acceleration 


z 

8 " “a v 

which is the initial acceleration response if unsteady-lift effects are 
ignored, and p is a mass parameter defined by 


* tm a 

An examinat ion of equation (25) indicates that for large values of 
the mass parmeter the alleviation factor K is substantially equal to 
the lift-response function k(s). Thus, any c han ges in k(s) are then 
directly reflected in similar changes in the normal-acceleration response. 

For smaller values of the mass parameter the integral equation has 
to be solved for K(s) by using operational methods (if k(s) and 
k^(s) are, or can approximately be, expressed by simple functions), by 
Iteration, or by numerical methods, numerical calculations have been 
made for p = 50 and 200 by using the two-dimensional unsteady- lift 
functions for incompressible and supersonic flow. The value of p = 50 
represents a tr ans port airplane flying at low or moderately high alti- 
tudes •, the other value represents a fighter or bomber airpla n e at low 
altitudes, or the tr ans port airplane at high altitudes. The results are 
shown in figures 19 to 22. 
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In using response functions for two-dimensional, flow to represent .. 
those for an actual tapered wing, basically two approximations are made, 
inasmuch as both the steady-state lift and the mann er in which it is 
approached as a function of time differ for a_ wing of finite and infinite 
span. The differences in steady-state values can be removed by usin g the 
unsteady-lift functions k(s) which are normalized with the steady-state 
lift; the differences in the manner in which this value is approached can 
be minimized by referring s to a more representative chord than the 
root chord, such as the mean aerodynamic chord. If s is redefined in 
this manner for the purpose of using the two-dimensional response func- 
tion to better advantage, the mass parameter p must also be referred 
to that chord rather than to c. Hie values of jj. and s used in 
figures 19 to 22 should, therefore, be considered to be defined in this 
modified manner. 

Comparison of the results presented in these figures with those for 
the unsteady-lift functions themselves indicates that, as expected, the - 
acceleration response follows the lift response not only for large values 
of p (around 200) but also during the first .few instants of travel (in 
which time the peak in the lift response due. to noncirculatory flow 
effects occurs in subsonic flow) for smaller mass ratios (around 50); 
however, thereafter the acceleration response does depend on the mass 
parameter. These statements are true for gusts both approaching the 
wing and overtaking it. For supersonic speeds calculations have been 
made only for a mass ratio p of 50, because even for this value the 
acceleration responses follow substantially the lift-response character- 
istics; they can be expected to do so to an even greater extent for 
larger values of the mass ratio. The conclusion can, therefore, be drawn 
from these curves that the acceleration response depends to a large 
extent on the magnitude and direction of the gust Bpeed. 

In the preceding paragraphs only sharp-edged gusts have been con- 
sidered. For traveling gusts, the intensity of which is a function of 
the distance behind the wave front d, the acceleration response can be 
determined by superposition from the sharp-edged gust response obtained ' 
by solving equation (25). Inasmuch as 



the gust intensity can also be expressed as a function of s, namely, 
W(s). In terms of this function and the acceleration-response function 
K(s), the normal acceleration for a time-dependent gust can be written 

as 
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( 26 ) 


The convolution process Indicated in this equation may tend to 
obscure some of the minor effects resulting from changes in gust speed, 
but the major conclusions- reached concerning the effect of such changes 
on the acceleration response to sharp-edged gusts are likely to be v ali d 
for the acceleration response to other types of gusts as well. 


Z(s) 


qSC. 


i' 


r 8 dwfs, ) 

W(o)K(s) + / K(s-S]_) — dB! 

J 0 ds 1 


RELATION BETWEEN BLAST WAVES AND TRAVELING GUSTS 


When a blast wave strikes a stationary object, the ins tantaneous 
pressures on the surface of. the object are proportional to the peak over- 
pressure of the blast, the constant of proportionality being a reflec- 
tion factor which depends on the geometric characteristics of the object 
and, for strong shocks, on the shock strength as well. (See ref. 13 .) 

This overpressure is due not only to the fact that the object 1 b initially 
exposed to the shock but also to the fact that it arrests the propaga- 
tion of that shock or deflects it. The manner in which the problem of 
calculating this overpressure (the diffraction problem) can be solved 
is indicated in reference l4, where the results of several such calcu- 
lations are also given. 

When the object is a wing flying initially through still air, Its_ 
response to the overpressure of the shock is similar, but in addition to 
this effect it also responds to the velocity behind the shock. The 
effects of overpressure and velocity overlap and c ann ot be divorced from 
each other readily. For instance, for a weak shock, which propagates 
substantially at the speed of sound, the overpressure p^ and, the 
velocity Wjj are related by the limiting form, for small overpressures, 
of the Rankine -Hugonio t equation, namely, 

Ph = paW b 

If this shock strikes a plate parallel to Its front, the initial pres- 
sure on its exposed surface Is 2^, whether the plate is stationary 
or moving parallel to the shock front, and the initial pressure coeffi- 
cient for the surface Is 
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However, this Initial pressure coefficient is also precisely the one 
associated with instantaneous entry into a vertical gust of intensity 
W-^; thus, consideration of either the overpressure or velocity yields 

a result for the initial response which includes the effects of the other 
in this case. On the other hand, after some time has elapsed, the over- 
pressure equalizes around the plate and, hence, produces no pressure 
difference directly, although in subsonic flow it can still influence 
the lift through the yorticlty shed while it was acting; thus, except 
for this induction effect in subsonic flow, the pressure difference is 
then due to the velocity effect alone. 

In this section only the effects of the change in relative velocity 
occasioned by entry into a weak blast are considered. In view of the 
preceding argument, these effects include at least some of the effects 
of the overpressure associated with the given blast. Inasmuch as atten- 
tion is confined to weak blasts, the effects of change in temperature 
(and hence speed of sound) and density across the blast wave, which pro- 
duce only second-order effects on the lift of the wing, will be disregarded. 

The relation between a blast wave and a traveling gust is indicated 
in the following sketches: 



The equivalent vertical gust Intensity is determined by a consider- 
ation of the conditions prior to and subsequent to gust and blast entry. 
Thus, 


W = Wfc sin 0 
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The other (horizontal) component of the hlast velocity is parallel to 
the chord and is equivalent to a horizontal gust which, within the scope 
of a linearized treatment, produces the same lift as a vertical, gust of 
intensity equal to 2ccq times the intensity of the horizontal gust. 
Thus, if this effect is taken into account. 


W = W- b (sin 0 + 2oq cos 0 ) 


(27) 


For "blasts striking from, almost straight ahead or behind, the second term 
may be more important than the first. 


The equivalent speed of the traveling gust is determined by a con- 
sideration of the time required for the gust and, the blast wave to trav- 
erse a given horizontal distance. Thus, the time required for the blast 
wave to pass over a given distance Ax along the flight path is 
Ax cos e/V b , and the time for the gust is Ax/V g . Therefore, the equiva- 
lent gust speed is 




8 cos 0 


V„ « 


8 cos 0 


( 28 ) 


inasmuch as for a weak blast wave, V-^ is approximately the speed of 
sound a. The equivalent value of V S associated with a blast must, 
therefore, always be supersonic. 


With the equivalent relations given by equations (27) and (28), the 
gust-response functions for supersonic values of V g arid acceleration- 

response calculations presented here can thus be used for the blast- 
penetration problem. The variation of the blast velocity W^ behind 

the blast front can be taken into account by expressing this variation 
as a function of b by using the relation 



cos 0 +■ 


and then by using equation (26) to calculate the acceleration response. 
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DISCUSSION 


The lift and. moment responses have heeh calculated herein for two- 
dimensional wings in incompressible and supersonic flow, for wide delta 
and rectangular wings in supersonic flow, and for very narrow delta 
wings. The reason for selecting these cases was that for them the func- 
tions of interest could be readily calculated. For three-dimensional .. 
incompressible flow and for two-dimensional compressible flow the func- 
tions k^_ and have been calculated in. previous investigations, but 
the methods used do not lend themselves to the calculation of the more . 
general response functions of interest here; also no simple methods are 
available for calculating the desired responses of wings with subsonic... 
leading edges in supersonic flow. 

The results of calculations for very narrow delta wings have been 
included despite the fact that the results of linear theory for these 
wings are of limited practical utility, I nasmu ch as for even relatively 
low angles of attack the wake of these wings tends to curl up and intro- 
duce deviations from linearity in force and moment responses. The 
reasons for including these results is that they are obtained very easily 
for incompressible' flow and serve as the only indication of aspect-ratio 
effects in subsonic flow available at present. In other words, their 
sign i ficance stems primarily from the fact that these wings serve as a 
theoretical limiting case. However, no numerical calculations have been 
made for these wings in compressible flow because such calculations 
would have required more effort than the results are considered to 
warrant . 

In this paper only entry into traveling gusts and blasts is con- 
sidered. However, the results also apply to . the problem of exit from 
traveling gusts a n d blasts, inasmuch as exit .from a gust of given 
intensity can be considered to represent an entry into a gust of equal 
and opposite intensity,* thus, this problem heed not be treated explicitly, 
and all statements made in this paper concerning responses due to entry 
into gusts and blasts are equally valid for exit from gusts anti, blasts . 

A full discussion of the implications and limitations of linear 
theory are beyond the scope of this paper. TEt should be pointed out that, 
in order for the results of this theory to be valid, all disturbances 
must be small; therefore, the gust intensity and the magnitude of the 
resulting motion must not be too large. For - most problems other t han 
those related to. blast loads, these conditions are likely to be satisfied, 
except possibly for very narrow delta wings, ' On the other hand, for 
blasts they may not be satisfied unless the airplane is fairly far away 
from the center of the blast. The limiting distance depends on the magni- 
tude of the blast; the orientation relative to the airplane (blasts from 
almost directly ahead of or behind the airplane result in low induced 
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angles of attack); the speed of the airplane (the higher the speed the 
smaller the angle.-of-attack change is caused by a given blast velocity; 
but for high subsonic, transonic, and high supersonic speeds, only s m a ll 
angles of attack can gener al 1 y be tolerated within linear theory); and 
the plan form (the aerodynamic forces on low-aspect-ratlo wings tend to 
be linear only in a small angle-of -attack range). Nonetheless, despite 
these limitations and despite the possible presence of the aforementioned 
pressure effects which have been ignored herein, the results presented 
herein which pertain to supersonic speeds of gust propagation should 
furnish a first step in the theoretical analysis of the blast-load 
problem. 


CONCLUDING REMARKS 


The lift and moment responses to penetration of sharp-edged traveling 
gusts have been calculated for wings in incompressible and supersonic 
two-dimensional flow, for wide delta and, rectangular wings in supersonic 
flow, and for very narrow delta wings. By using the two-dimensional 
indicial-lift functions, some calculations of normal-acceleration response 
have been made for two mass ratios. 

The results of these calculations indicate that the propagation 
speed of the gusts has a large effect on the lift- a n d moment-response 
functions. For incompressible flow (and presumably for compressible 
subsonic flow as well), peaks exist In the early portion of the lift 
response which may be much larger than the steady-state value. Some 
peaks also occur in the lift-response functions for supersonic speeds 
but are much less pronounced and exceed the steady- state value in only 
a few instances. 

The normal-acceleration responses tend to follow the lift response 
very closely In the first few Instants of penetration provided the masB 
of the airplane is not impracticably small; thus, the large peaks which 
exist In the lift response at subsonic speeds are duplicated In the 
acceleration response. After the first few instants the acceleration 
response depends to a relatively larger extent on the mass of the air- 
plane. 

The relation between gusts traveling at supersonic speeds and blast 
waves has been indicated, and the manner in which the calculated lift and 
moment responses can be used in a linearized approach to the blast-load 
problem has been outlined. 


Langley Aeronautical laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , January l6, 1957* 
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APPENDIX 


RELATION BETWEEN UNSTEADY -LIFT FUNCTIONS FOR GUST PENETRATION 
AND FOR INDICIAL FLAP DEFIEOTION 


As pointed out in the text of this paper, the unsteady-lift or 
unsteady -manent functions for gust penetratiog. cannot, in general, he 
obtained directly from the lift and moment functions for lndicial angle- 
of -attack change. However, they can he obtained from the unsteady-lift 
and unsteady-moment functions for indicial flap deflection. This approach 
is similar to the one used by Ktissner in reference 1 for stationary gusts. 
However, in reference 1 consideration is given initially to sinusoidal 
gusts and flap deflections; whereas here the. desired results are estab-.. 
lished by working with indicial responses directly. 

As may be seen from equation (3), the boundary condition for gust 
penetration is, at any instant t', the same as that for a wing with a 
leading-edge flap which extends to the point x' = (v + Vg)t'. Conse- 
quently, if the lift or moment due to indicial flap deflection is known, 
the lift due to gust penetration can be obtained by superposition. Thus, 
if the normalized unsteady-lift function for Indicial deflection of a 

trailing -edge flap extending rearward from x' «= g ^ is designated by 

kg(s,g), "the unsteady- lift function for gust penetration can be obtained 

by considering segments of the wing to deflect indicially in a progres- 
sive fashion; thus, 

k(s) « [lsg(s,0) - kg(s,A|)] + [kB(snAs,A£) - k 6 (s-As,2Ag)] + • ■ • 

» -Ag (kg(s,0) + kg(snAs,Ag) + kg(s-2As,2Ag) + . . .] 
with — _ — 


Ag 


V 4- V, 


£ As 


V 


_ As 
* 
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Hence, in the limit as A| -»0, 



This relation pertains to stationary gusts (A = l) as well as for 
all gusts traversing the wing from the leading edge to the trailing edge 
(A > -l). For gusts traversing the wing from the trailing edge to the 
leading edge < -l), the following relation may be obtained in the 
same manner; 



Equations (Al) and (A2) represent the deBired relations between the 
unsteady-lift functions for gust penetration and indicial flap deflection. 
Identical relations exist between the unsteady-moment functions. 

The unsteady lift (or moment) due to an indicial flap deflecting, 
as used in the preceding relations, is the lift on a wing which has a 
flap deflected by a unit angle, is initially at rest, is accelerated 
instantaneously to speed V at time t' = 0, and then continues to fly 
at that speed with zero angle of attack. This lift is related to but 
not necessarily identical to the lift on a wing initially flying at 
speed V an d experiencing an abrupt flap deflection at time t* = 0, 
because such a deflection involves transient effects, namely, those due 
to an impulsive rate of rotation (a condition which implies an impulsive 
vertical velocity proportional to the distance from the axis of rotation), 
which are spurious for the present purpose. 

Similarly, the lift for a flap which is abruptly given a constant 
rate of vertical displacement at time t’ = 0 contains other effects 
which are spurious for the present purpose, namely, those due to the 
chordwise discontinuity in displacement at the leading edge of the flap 
which must exist in this case. 

The unsteady- lift (or moment) functions for abrupt flap displace- 
ments can be obtained from those due to oscillating displacements, which 
are often calculated for use in flutter studies, by means of the relation 
(see ref. 5) 



(A3) 
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where Fg is the real part of the normalized unsteady-lift function for 
sinusoidal flap displacements, and k is a reduced frequency equal to 
A similar relation exists between the moments . 

The following two examples indicate the manner in which the required 
functions may be obtained. 

For supersonic two-dimensional flow, a trai ling-edge flap of chord 
Cf experiences the same lift and, moment as If it were a wing with chord 

Cf j and the rest of the wing contributes nothing to the lift and moment. 
Hence, 

kg(s,|) = ^ bJ 

where k]_(s) represents the unsteady-lift function for indicial angle- 
of -attack change, and inasmuch as £ = 2^1 - , 

* 6 ( 8 , 6 ) = |(2 - (&) 

Similarly, if C^iXs) is the unsteady -moment function for indicial 
angle -of -attack change and-the moments are taken about the leading edge_^ 

<^(=,6) - i(2 - - |(2 (A5) 

where C]^, the lift-curve slope, is the value with which k^(s) and 
k(s) are normalized. 

For incompressible two-dimensional flow the normalized complex lift 
function for a tralling-edge flap undergoing oscillatory displacements of 

amplitude Hq£ with frequency cn can be obtained from those of the terms 

of equation (22) of reference 15 which contain the product = H. The 

following result is obtained in this manner: 



2 c H 


2 V 



e H 
2 VJ 


where C is a function of the reduced frequency k which is defined in 
reference . 15 > and generally is known as the jjheodorsen function, and where 
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the functions 3^, T^, and Tgi are defined in reference 15 in terms- 

of a parameter c which in the notation of the present paper is | - lj 
thus, 

% - ( 8 - 1)^2£ - I 2 - cos _1 ( 1-1) j 

\ (A 6 ) 

TqO = ^2| - I 2 + cos-M | - 1) 


In the expression for C L the terms in the first bracket are asso- 
ciated with noncirculatory flow, the others with circulatory flow. In 
each bracket the first term is due to the vertical displacement of the 
flapj therefore, it is not pertinent to the present purpose. (Only the 
contribution of the potential given in eq. ( 7 ) of ref. 15 is pertinent 

here; the potential given by eq.. ( 6 ) of ref. 15 ia not pertinent.) 

Hie desired normalized steady-lift function is, therefore, for unit down- 
wash, that is, for 3s. = - 1 , where 

%(k, 8 ) = - + C (k) ^ (AT) 

Thus, equation (A3) yields 

%U,|) -§t 8 (.) CAB) 

where 8 (s) is the unit impulse (Dirac delta) function and k-^ c (s) is 

the continuous part of the judicial- lift-response function, which is the 
part obtained directly from the circulatory part of the lift due to sinus- 
oidal vertical oscillations, namely, C(k), by using equation (A3). 

Similarly, the moment about the quarter-chord point can be obtained 
from equation ( 23 ) of reference 15 by setting a = - jy and retaining only 
the terms containing [31; therefore, 

tv - - + Tg 1 1 + |] 

This moment iB due entirely to potential effects. Again, the first term 
and part of the second are due to the chordwise discontinuity and are 
spurious. The other two terms can be expressed in terms of the functions 
1^ and T 10 used for the lift and a function Tg defined by 
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3 ^ 


T 8 = - |(2| - | 2 ) 3/2 - (£ - 1 )\ (A9) 

Thus, again for - -1, this moment function becomes 

C^(k, s) . - + I w ) + f lig + ( 5 - 1) J (AID) 

Hence, hy using equation (A3), the following expression is obtained for 

CmgCs,!): 

QtngC 0 , 2) - - |(^ + T 10 )i(s) +-i [t 8 + (5 - S(s) (All) 


The manner in which these functions can .he -used to calculate the 
lift and moment responses to traveling gusts by means of equations (Al) 
and (A2) is illustrated in the text of this paper for incompressible 
two-dimensional flow. -- 
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fcABL* 1.- OHfUADT-Lira TOJCHOB k(a) KB nOCMHBSS 
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H»-DDIHrB33QWLL TUN 


Guat condition 

Bangs of a fop - 

0<«S £|X| 

w 

■ * a|X| 

Quart upporoachlng vlng and gust 
recoding m vlng at spools laaa 
than that of vlng (X > 0) 

+ j^aoe -1 * 8 - + \/s(SX - a) - 

B / X + ■ 1 a(2X - a) + MX - a) 

M\^ - X - a| ^ XX 


Oust overtaking vlng (X < 0) 

i^L + ^ooa“ 1 ^-±-i+i^l/|a(BX+ a)| - 
2 ^Ib j X + 2X1 otn -l »(2X + ■) + X(X + a) 

1 , ^ 1 1 1 a) 


*In the hurt teas the function oob"^- 1 a raplaood by ooal>~^ vban a % 2|X| - V. 


*abi* a.- TBBnunz-janR pokbcch c^(b) rat jkxnebbbxbii nM-minatau, tun 

Qtcnsnt la taken about quarter -chord point] 


Goat ifnH-MrM 

Bangs of a for - ^ 

0 < a S 2|X| 

■ i 2 |x| 

Oust approaching vlng or reoedlng 
fraa vlng at apooda laaa than 
that of wing (X > 0) 

■ 

(S 

*4 

K 

1 

ri 

ft 

1 

0 

Guat overtaking vlng (X < 0} 

-fa^A)(gx + s)\fs|ex + a| 
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SABI* 3— IHBS8CB* C UUnCX tt t Cp TOE ajRHJBDUC WO-DMOBIOIAL TIM MB) IflMUBG 
C UBff ili-Lim l/q FOB OEUtA VJXB KCffl 0UPBJ9OCC ISADXR FUSS 

(s) Ousts app r oa ch i n g ving ar receding trcm. sing at speeds 
Isas than that of ving (Mg > -M) 


Ooeffioiont t - Oj 
X«t| 
x - c - Mt 


Begton In x,t plane for - 


x ■ tj x - 


1 Mr + t 


Mr 


t 6 PTI 


* 

ate + M(x 4- Mt) coB . 1 Mr + t + 

Mg 2 _ l -jtP - 1 x + Mt 

T — } ■ - — [x(nJ + M) + MgtflL 2 )! - Mg - »)] cos -1 * ~V . 


Wm 2 - 1 




In the last t«m the cos -1 is replaced by cosh -1 for | Mg] < 1. 

i 

’■^srrr#^) (*>V.-»<V< -*) 

5 - + K ) + v k* - "H - *8 (-<VJ«8>3) 


i.a&H fr + ift ) + 

i VM .r^ — r 


x(mb 3 + m) + Mgtfoftt - ife - at)' 

L.a 3/2 


i (x :> Mgt; J< < Mg < -l} 
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(a) Ousts approaching wing or receding from vlng at speeds 
below that of wing (Mg > -M) 
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TABIE 7.- UNSTEADY-LIFT FUNCTION k(a) KE A NAEROW-DKITA WING IN COMPKESSIBIB FLOW 
(a.) Ousts approaching wing or receding from wing at speeds lass than that of wing 


of s far - 


Mg > -H 
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(h) Gusts overtaking wing [Mg < -It) 
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xKADx PTKIBXHQ tCMBT ABOTO TBs LSS331W EDGE TOR A HAHFDW QEUCA WIBG IE CCMFI 
(a) Gusts ajrproeabing wing ar recoding from vlng st speeds Uses .than that of wing 
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(b) Gusts approaching from tr ailing edge. 
Figure 1.- Concluded. 
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Slrtanc* MnlM, •> MMlctaorda 


(a) Gusts approaching from leading edge. 

Figure 2 .- Indicia! -nranent functions for a wing in incompressible two-dimensional flow when 
penetrating trave ling sharp-edged gusts (moment about quarter chord) . 
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(a) Gusts approaching wing or receding from wing at speeds less than 

that of wing. M g > -M. 



Figure 3 .- Sketches of geometric characteristics for wings in two- 
dimensional supersonic flow and for delta wings with supersonic 
leading edges. 
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Dlatanca traveled, a, awl chorda 

(h) Gusts approaching from trailing edge 

1 Figure 4.- Concluded. 
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Dlstanca trwalad, a, root — lclinrda 


(a) Gust approaching from leading edge. 

Figure 7 .- Indicia! -moment functions for a wide delta wing flying at M *» 2 

traveling sharp-edged gusts (moment shout leading edge) 


when penetrating 







Dictum trml*4, I, *mlt8iot6a 


(a) A = 1. 

Figure 9.- Indicial-lift functions far vide rectangular vlngs of various aspect ratios when 
flying at M = 2 and penetrating traveling Bhaj^p -edged gusts approaching from lead i ng edge 
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(b) A - 2. 

On 

Figure 11.- Continued. -J 
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Distance t rare led, s, antrl chorda 


(a) A - 1. 

Figure 15 .- Indlelal -moment functions for vide rectangular wings of various aspect ratios when 
flying at M = 2 and penetrating traveling' sharp-edged gusts approaching from leading edge 
( m oment about leading edge) . 
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ntatauc* tr***l*d, a, sadcborda 


(a) A = 1. 

Figure 15.- Tnfli - mrmpr't . functions for wide rectangular wings of various aspect ratios when 
flying at M =» 2 penetrating traveling sharp -edged gusts approaching from trailing edge 

(moment about leading edge) . 
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Distance traveled, ■, root semlchordi 


(a) Gusts approaching from leading edge. 

Figure 17 •- Indicial-lift functions for very narrow delta wing in 
incompressible flow when penetrating traveling sharp -edged gusts. 
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Dlrtanca tmalad, *, root —lehords 


(a) Gusts approaching from leading edge. 


Figure 18. - Indicial-mcanent functions for very narrow delta wing In incompressible flow when 
penetrating traveling sharp-edged gusts (moment shout two-thirds chord) . 
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Distance traveled, a, sealcbords 

Figure 19 .- Alleviation factor for a wing "with mass -density ratio of 
H = 50 in incompressible two-dimensional flow when penetrating 
traveling sharp -edged gusts approaching from leading edge. 
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Figure 21.- Alleviation factors for a wing’ with mass -density ratios of 
= 50 and 200 in incompressible two -dimensional flow when pene- 
trating a traveling sharp -edged gust approaching from trailing edge 
at = — 0 . 23 • 
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Slrtuce trav»X= 


Figure 22.- Alleviation factor for two-d 
with a mass -density ratio of |i = 50 
sharp -edged gusts approaching from. ] 
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